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Abstract 

In a 3-wave mixing process occurring in some nonlinear optical medium when ajmode inter- 
acts with both a3,mode and ajmode, we theoretically study the squeezing effect generated by the 
operator S3 = exp[/x(aia2 — a\a\) + v{a±az — ajaj)]. The new 3-mode squeezed vacuum state in 
Fock space is derived, and the uncertainty relation for it is demonstrated, It turns out that S^may 
exhibit enhanced squeezing. By virtue of the technique of integration within an ordered product 
(IWOP) of operators, we also derive S3 's normally ordered expansion. The Wigner function of new 
3-mode squeezed vacuum state is calculated by using the Weyl ordering invariance under similar 
transformations. 

PACS 42.50.-p - Quantum optics 

PACS 03.65.-w- Quantum mechanics 



1 Introduction 

Nowadays quantum entanglement is the focus of quantum information research and attracts many 
interests due to its wide applications in quantum communication 0] [2] . Entangled states have 
brought much attention and interests of physicists J31 |4j . The usual two-mode squeezed state, gen- 
erated from a parametric amplifier flSJ, not only exhibits squeezing, but also quantum entanglement 
between the idle-mode and the signal-mode in frequency domain. Therefore, it is simultaneously a 
typical entangled state of continuous variable. Theoretically, the two-mode squeezed state is con- 
structed by acting a two-mode squeezing operator S 2 = exp[A(aia2 — <4c4)] E[3 on the two-mode 
vacuum state |00), i.e. S 2 1 00) =sechAexp[— a\a\ tanh A] 1 00), where A is a squeezing parameter, and 
Oj(a]) Bose annihilation (creation) operator satisfying [04, oj] = <5y. Using the relation between Bose 
operators (oj, aj) and the coordinate and momentum operators 

one can recast S% into the form 

S 2 =cx P [iA(Q 1 P 2 +Q 2 P 1 )], (2) 

noting 

[QiP 2 ,Q 2 Pi] = HQ2P2-Q1P1), 

[QiP 2 ,±{Q2P2-QiPi)] = 2Q 1 P 2 , (3) 
[Q2P1, i (Q2P2 - QiPi)\ = -2Q 2 P U 
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thus there involves a SU(l, 1) algebraic stricture. In the state 52 |00), the variances of the two-mode 
quadrature operators of light field, 

+ (4) 

satisfying the commutation relation [X, — |, exhibiting the standard squeezing, i.e., 

(00| S 2 f X 2 S 2 |00) = ie- 2 \ (00| Sly 2 S 2 |00) = ie 2 \ (5) 

which satisfy (A£)(A«p) = i 

An interesting question naturally arises: if o| mode in a nonlinear optical medium, interacting 
with both a\ mode and a\ mode (e.g., a three-wave mixing), and the corresponding three-mode 
exponential operator is introduced as 

S3 = exp[/i(aia 2 — a\a\) + ^(aia 3 — <4t4)]. (6) 

Using Eq. (Q]) we can recast S3 into the form 

S 3 = exp [i/i (Q 2 Pi + QxP 2 ) + Iv (Q3P1 + Q1P5)] , (7) 

where /x, are two different interaction parameters, then what is its squeezing effect for the 3-mode 
quadratures of light field? 

To answer this question we must know what is the state S3 |000) (|000) is the 3-mode vacuum 
state) in Fock space, for this aim, we should know what is the normally ordered expansion of S3. 
But how to disentangle the exponential operator S3? Because there is no simple SU(1, 1) algebraic 
structure among Q2P1, Q1-P2, Q3P1 and Q1P3, the disentangling seems hard. Thus we turn to appeal 
to Dirac's coordinate representation and the technique of integration within an ordered product 
(IWOP) of operators Il8l [9l [lOl I11H to solve this problem. Our work is arranged as follows: firstly 
we derive the explicit form of S3 |000) , then we demonstrate that it really satisfies the Heisenberg 
uncertainty relation and may exhibit squeezing enhancement. We also employ the technique of 
integration within an ordered product (IWOP) of operators to derive the normally ordered expansion 
of S3. The Wigner function of S3 |000) is calculated by using the Weyl ordering invariance under 
similar transformations H121 [T3l [l4ll . 

2 New 3-mode squeezed vacuum state 

For the sake of convenience, we rewrite S3 in Eq. (0 as the following compact form, 

S 3 - oxp i^A, ,/>, .'•../ = 1, 2, 3, (8) 
where the repeated indices imply the Einstein summation notation, and 



A = m , (9) 




thus 



its inverse is 



cosh r cos 9 sinh r sin 8 sinh r 

+ cos 2 9 cosh r ^ 
i29 1 „t, „ 1 \ .2 , 



A -' cosflsinhr sin 2 9 + cos 2 9 cosh r (coshr - 1) ), (10) 



sin 9 sinh r (coshr -1) sin 6> coshr 



cosh r — cos 9 sinh r — sin 9 sinh r 

e ~ A = I — cos 9 sinh r sin 2 9 + cos 2 9 cosh r sm29(cosh? -1) j q-q 

- Sin 9 Sinh r sm2g(coshr-l) g - n 2 q QQgJj r _|_ cog 2 Q 



where we have set 

M . a _ v 



y/^T^,cos9= sin 6= -, (12) 



r 



noting that A is a symmetric matrix. Using the Baker-Hausdorff formula, 

e A Be- A = B+[A,B} + ±[A, [A, B]] 

+ i[A,[A,[A,B}]} + -- - , (13) 

we see that 53 causes the following transformation 

S^QkSs = (e- A ) kl Q u S^PkSs = (e A ) H P t . (14) 
It then follows that S 3 ~ 1 akS3 = {e~ XK )uiCLi, i.e. 

i 4. 4. 

S 3 ai53 = a\ cosh r — a 2 cos sinh r — a 3 sin sinh r, 
S^~ 1 a2<S , 3 = — a{ cos sinh r + a2 (sin 2 9 + cos 2 cosh r) 

+ ia 3 (cosh r - 1) sin 20, (15) 

t 1 
5*3 &3S3 = —a} sin sinh r + — a 2 (coshr — 1) sin20 

+ a 3 (sin 2 9 cosh r + cos 2 0) . 

Noticing that S3 = S^ 1 and S3 (n, z/) = .S3 (— ju, —v), from Eq. (fT5l) we also have 

6*301 S3 1 = ai cosh r + a\ cos 9 sinh r + 4 sin sinh r, 
S3C12S3 1 = a\ cos0sinhr + a 2 (sin 2 9 + cos 2 coshr) 

+ ~a 3 (coshr - l)sin20 (16) 

S3CL3S3 1 = a\ sin 9 sinh r + —a 2 (coshr — 1) sin 29 
+ 03 (sin 2 9 cosh r + cos 2 9) . 

For convenience to write, we set S 3 |000) = 1 1 000) . In order to obtain the explicit form of ||000), 
using Eq. (ll5D and a\ |000) = 0, we operate a\ on 1 1 000) and obtain 

a x ||000) = SaS^aiSs |000) 

= S3 (ai coshr — a\ cos 9 sinh r — a\ sin sinh r) 1 000) 
= 5 3 (-4 cos0sinhr- a\ sin sinh r)^ 1 ^ |000) 

= -S* 3 (4cos6lsinhr + 4sin6»sinhr)5 3 " 1 ||000) , (17) 
then we continue to use Eq. dl6D to derive 

ai || 000) = — {[a-i cos 9 sinh r + a 2 (sin 2 9 + cos 2 9 coshr) 
+ — (cosh r — 1) sin 29] cos 9 sinh r 

+ [ai sin sinh r + —a 2 (cosh r — 1) sin 20 

+ aj (sin 2 cosh r + cos 2 0) ] sin 9 sinh r} || 000) 
= — (ai sinh 2 r + a\ cos cosh r sinh r 

+ i4 sin6,sinh2r )ll° 00 )' ( 18 ) 



so we reach the equation 

||000) = -tanhr(4cos6» + 4sin6»)||000). (19) 
Similarly, operating a 2 on 1 1 000) and using Eqs.( 115D and ( |16t yields 

a 2 ||000> = SaS^azSa |000> = S 3 (-a\ cos6»sinhr) |000) 
= S 3 (-a\ cos (9 sinh r)^ 1 ||000) 
= — (a\ cosh r + 0,2 cos 9 sinh r 

+ a 3 sin 9 sinh r ) cos 9 sinh r \ | 000) , (20) 

which leads to 

[a 2 (1 + cos 2 6» sinh 2 r) + ^a 3 sin 26> sinh 2 r] 1 1 000) 

= --a| cos 6> sinh 2r 1 1 000) . (21) 
On the other hand, operating 03 on 1 1 000) and using Eqs. (115D and (1161) yields 

a 3 || 000) = S 3 5^ 1 o 3 S r 3 |000) = 5 3 (-aJ sin6»sinhr) |000) 
= ^(-alsinflsinhr)^ 1 ||000) 
= — (a| cosh r + 02 cos 9 sinh r 

+ a 3 sin 9 sinh r ) sin 9 sinh r 1 1 000) , (22) 

i.e., 

[03 (l + sin 2 9 sinh 2 r) + ^a 2 sin 26* sinh 2 r] || 000) 

= -~a\ sin 6» sinh 2r 1 1 000) . (23) 

Combining Eqs. (l2lT) and (|231) we have 

a 2 ||000) = -aJtanhrcosflHOOO) , (24) 

and 

a 3 || 000) = -oj tanhr sin (9 ||000) . (25) 
From Eqs.(fT9l).(|24l) and (|25|) . we may predict that 1 1 000) has the following explicit form: 

1 1 000) = N exp [- (4 cos 9 + o| sin (9)4 tanh r] 1 000) , (26) 

where N is the normalization constant, which can be determined by (000 1 1 000) = 1, and we 
calculate N — sechr. 

3 Squeezing property and quantum fluctuation in 1 1 000) 



Squeezing is an important phenomenon in quantum theory and has many applications in various 
areas in quantum optics and quantum information II 1 5 II . In this section, we examine the quadrature 
squeezing effects of || 000) . The quadratures in the 3-mode case are defined as 



which satisfy the relation [X%, X 2 ] = §. Their variances are (AXi) 2 = (X?) — (Xi) 2 , i = 1,2. Noting 
the expectation values of X x and X 2 in the state || 000) is (Xx) = (X 2 ) =0. With the help of Eq.([l~5l). 
we can calculate that the corresponding variances in the state 1 1 000) : (noting A is symmetric) 

(AXi) 2 = (000 1 S^XfSa 1 000) 

1 3 3 

= ~ E E( e ~ A ) fc *( e ~ A )^ < 00 °i Q^ 1 

i=l 3=1 
3 3 

= ^ E X>~ A M e ~ A );< 

»=i i=i 

-^EE( e " A )^( e " A ^ 5 « 
»=i j=i 

^E(e- 2A ),3, (28) 



12 



and 



(AI 2 ) 2 = (000| 5 3 - 1 X 2 5 3 |000) = ^E(e 2i %- C29) 

The explicit form of the matrices e 2A and e~ 2A can be derived from Eq. (1101) and dill) , so we can 
obtain 

9 1 

(AX t ) = — [(2cosh2r + 1) + sin 26>(cosh 2r - 1) 

+ 2 (cos 9 + sin 0) sinh 2r] , (30) 



and 



9 1 

(AX 2 ) = — [(2cosh2r + 1) +sin26>(cosh2r - 1) 

- 2 (cos 9 + sin 0) sinh 2r] . (31) 



We can successfully verify 

(Ali)(AJ 2 ) 



= V ( 4 cosh 2r + 4) + (l - 2 sinh 2 r sin 26 
^ ^Y // ( 4cosh2r + 4 ) + (l-2sinh 2 r) 2 



= — \/-cosh4r + — ^ -, (32) 
12 V 2 2 4' 

which confirms the uncertainty relation of quantum mechanics. 

To see the trend of squeezing effects in the X x — or ^—direction, we plot (AX^ 2 and (AX 2 ) 
as the function of parameter fi for different v in Fig.l. When v = 0, it exhibits the usual two mode 
squeezing effect depending on the varying fi, (AXi) 2 increases accompanying (AX 2 ) 2 decreases; 
when v = 0.5, (AXi) 2 increases more than the case of v = 0, which exhibits enhanced squeezing 
in certain domain of fi. In Fig.2, we plot the uncertainty value (AXi)(AJf 2 ) as the function of r for 
different 9. 



4 Normally ordered form of S3 

We calculate the normally ordered form (denoted by : : ) of S3 by inserting the completeness rela- 
tion of coherent state 

f d 2 Zl d 2 z 2 d 2 z 3 

S3 = / 5 \z\z 2 z 3 ) (ziz 2 z 3 \ . (33) 



where \ziz 2 z 3 ) is the three-mode coherent state and \zi) = exp[— — h Zio\] \0i), i = l,2, 3. 
Using the relations in Eqs.( 115D and ( 1161 ), we have the explicit relation of S3 \ziz 2 z 3 ) 



)*2 



S3 \ziz 2 z 3 ) 

3 1 ,2 

= exp(- -^~) S 3 ex P ( 2 i4 + z *4 + *3<4) S3 1 Si |000) 

i=i 

1 I Zi 1 2 

= : — exp(— } — h ziz 2 cos#tanhr + ziz 3 sin#tanhr) 

cosh r ' 2 

i— 1 

x exp{ [at z\ + (at (sin 2 9 cosh r + cos 2 9) at (cosh r — 1) sin 29); 

coshr 1 2 J 

+ (ag (sin 2 + cos 2 6* coshr) — —a\ (coshr — 1) sin2C)z3]} 

x exp[-a|tanhr(4cos6> + (4sin6l)] |000) . (34) 

Substituting Eq.d34"l) into Eq. (l33l) , noticing that |000) (000| = : exp(— a\ai — a\a 2 — a\a 3 ) : , and using 
the following formula 

-iexp(c|z| 2 +ez + ^*) =-^e-^, Re(C)<0, (35) 

as well as the IWOP technique, we can obtain the explicit normally ordered expansion of S 3 : 

S3 = expl— at (at cos 9 + at sin#) tanhrl 

cosh r J ' 1 

1 — cosh r t t 2 

x: exp (a \a\ + a 2 a 2 cos 9 

coshr 

+ 0303 sin 2 9 + — a 2 a 3 sin 29 + —a 2 a 3 sin 29)} : 

x exp[ai(a2 cos 9 + ai sm9) tanhr]. (36) 

5 Wigner function of 1 1 000) 

Wigner distribution function of quantum states H161 1171 [18ll is widely studied in quantum statistics 
and quantum optics and is very important tool for a global description of nonclassical effect in the 
quantum system, which can be measured by various means such as photon counting experiment and 
homodyne tomography. Now we derive the Wigner function of 1 1 000) by using a new method. 

Recalling that in Ref. [11411 we have introduced the Weyl ordering form of single-mode Wigner 
operator Ai (qi,pi), 

Ai(gi.pi) = : <5(<zi -Qi)S( Pl -P x ) \ (37) 



where the symbols ; ; denote the Weyl ordering, while its normal ordering form is 

1 



Ai (qi,Pi) = - : exp - (q 1 - Q^ 2 - (pj - P x f 

IT 



(38) 



Thus the Wigner function for 1 0) can be easily expressed as (0| Ai (<?i,Pi) |0) = - exp(— q\ — p\). Note 
that the order of Bose operators ai and a\ within a normally ordered product (or a Weyl ordered 
product) can be permuted. That is to say, even though [oi,a|] = 1, we can have : a\a\ : =: a\a\ : 
and;aia^ = '. a \ a i '.- The Weyl ordering of operators has a remarkable property, i.e., the Weyl- 
ordering invariance of operators under similar transformations, which means 



u\ (ooo) \u~ l = ^(ooo)r 1 ;, 

as if the "fence" ; ; did not exist when U operates. 

For 3-mode case, the Weyl ordering form of the Wigner operator is 



(39) 



A 3 (q,p) = '5(q-Q)5(p-P) 



(40) 



where Q = (Qi, Q 2 , <2 3 



(Pi,P2,P 3 ) T , q= (qi,q2,q3) T and p = (pi,p 2 ,p 3 ) T - Then according 



to the Weyl ordering invariance under similar transformations and using Eq.( [14t , we have 

S^A 3 (q,p)S 3 

= ]S(q k - (e- A ) fcl Q 4 ) 5 (p fc - (e A ) fe4 P) 
= [ 8 {{e h ) kl q, - Q k ) 5 {{e~%, Pl - P k ) 

= ;*( q '-Q)*(p'-p); 

= A 3 (q',p'), 

where q' k = (e A ) fei g i , p' k — (e~ A ) ki pi. Thus the Wigner function of || 000) is 

(OOOlSj'Aa (q,p)5 3 |000) 



(41) 



= (000|A 3 (q',p')|000) 



T„2A, 



7t exp I — q e q — p e 



,T„-2A 



p) 



(42) 



where e 2A and e 2A are given by e A in Eq. dlOD and e A in Eq. dllD . respectively. 

In summary we have shown that the operator S3 = exp[^t(aia2 — a\a\) + ^(aia 3 — aja^)] is a 
new 3-mode squeezed operator by calculating the quantum fluctuation for 3-mode quadratures. We 
have obtained the new 3-mode squeezed vacuum state and derived the normally ordered expansion 
of S3. The IWOP technique brings convenience in our derivation. 
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Figure 1: (Colour online) The quantity (AXi) 2 and (AX 2 ) 2 as the function of squeezing parameter 
H for different case v = and v = 0.5. 




Figure 2: (Colour online) The uncertainty value (AXi)(AX 2 ) as the function of r for different 0. 



